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Abstract 

We compute the effective potential Kfr(v^) for one-component real scalar field in 
three Euclidean dimensions (3D) in the case of spontaneously broken symmetry, from 
the Monte Carlo simulation of the 3D Ising model in external field at temperatures 
approaching the phase transition from below. We study probability distributions of 
the order parameter on the lattices from 30'^ to 74^, at L/^ ^ 10. We find that, in 
close analogy with the symmetric case, ip^ plays an important role: Kfr(v') is very 
well approximated by the sum of (p"^, cp^ and (p^ terms. An unexpected feature is the 
negative sign of the (p^ term. As close to the continuum limit as we can get ~ 7.2), 
we obtain 

Ces ~ ^d^^d^ip + l.7{ip^ - r]yiip' + ri^). 

We also compute several universal coupling constants and ratios, including the combi- 
nation of critical amplitudes C~(/f )~^i?~^. 
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Introduction 

This work continues our previous Monte Carlo study which was devoted to the 
effective potential Vcs{<f) for the theory of one-component real scalar field in three Eu- 
clidean dimensions in the symmetric (paramagnetic, PM) phase. There we have found 
that Ves is very well approximated by the sum of v?^, and (p^ terms, and computed 
universal 4-point and 6-point couplings. Here we turn to the spontaneously broken 
(ferromagnetic, FM) phase of the same theory. To compute the effective potential in 
the FM case we use largely the same approach we have developed for the PM case, 
with some modifications. As the detailed discussion, with all necessary references, was 
given in m , we describe it here only briefly, concentrating mostly on the points specific 
for the broken phase. 

Monte Carlo computation 

We consider the Ising model in external field, 

{M <ij> i 

on a simple cubic lattice of the size (from 30'^ to 74'^) with periodic boundary 
conditions. The critical coupling is [3c ~ 0.221655 0. We use Swendsen-Wang cluster 
Monte Carlo algorithm to generate the Boltzmann ensemble of configurations, for 
given coupling (3 and external field J. For every configuration we measure the order 
parameter (magnetization per site) 

1 ^ 

<^=T7E'^- (2) 

i=l 

where N = IJ^ = Vl is the total number of sites. Thus we obtain probability distribu- 
tions P{^p), in form of histograms (Fig. |1|). The relation between P{(p) and Ves{<p), for 
sufficiently large volume Q, is 

P(<^) oc 7Kff(<^)exp {-QV,M + QM . (3) 

Now one can check whether it is possible to fit the set of probability distributions for 
given P and several values of J with this formula, using this or that ansatz for Ves{(p). 

We see that the situation is analogous to that in the PM phase the ansatz 
Ves{(p) = np"^ + uip^ works poorly, while 

Kfr(<^) = V + V + (4) 

provides a perfect fit (Figs. |l], 0). In the remaining part of the paper we discuss 
extraction of universal quantities from r, u, w obtained from such fits. 
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Magnetization 

Figure 1: Probability distributions P{(p) for magnetization per lattice site, for the Ising 
model ([![). The points are from Monte Carlo simulation; the solid line is the fit with 
Vefi{ip) = rif^ + uip^ + wip^ (all 5 histograms are fitted simultaneously with one Kffj- 



60 




Magnetization 

Figure 2: The Monte Carlo points are the same as at Fig. |I|. The solid line is given by 
(Qj, Veff{(p) = rip"^ + u(p^, where r and u are fixed by the requirement that spontaneous 
magnetization and zero field susceptibility are correctly reproduced. 



2 



Problems with fitting procedures 

In the FM phase several comphcations arise, that make determination of parameters 
of the effective potential much more tricky than in the PM phase. 

First, it turns out that the ratio of lattice size to correlation length L/^ has to be 
much larger in the FM phase than in the PM phase, to get finite size effects under 
control. The formula @ is valid in the limit of large L/^. It turns out that to get 
finite volume corrections to as small as they are in the PM phase at L/i^ = 4, one 
has to go to L/,^ as large as 10 in the FM phase. 

Thus it is difficult to go to large S,-, as the lattice size becomes prohibitively large. 
This, in its turn, makes it very difficult to study the finite cutoff effects and to extrap- 
olate to ^ — >■ oo. While it was possible to do such an extrapolation, with a reasonable 
degree of confidence, in the PM case, where we could study ^ as large as 14.1, in the 
FM case we can only get a qualitative estimate of finite cutoff corrections, and have no 
choice other than to take the values from our largest lattice, without attempting any 
extrapolation. 

Secondly, the applicability of eq. (|^) for various values of ip depends on the ratio 
L/^{ip), where ^{(p) is determined by V^^{(p), and deteriorates quickly at \ip\ < M, 
where M is spontaneous magnetization. Thus we have to use for the fitting only the 
subset of data corresponding to \ip\ > M, losing considerable amount of information 
(practically, a half of the J = histogram). 

So one has to take into account that 

1) Including larger values of J into the fit reduces statistical errors, but increases 
systematic errors connected with the finite cutoff (i.e. the smallness of ^), as larger J 
mean smaller ^. 

2) It would be nice to use the information contained in the whole J = histogram, 
including \ip\ < M. 

One can see how this works, from Table |1|, where we have compiled some observ- 
ables that can be measured directly, as well as position M of the minimum of Ves 
(that corresponds to spontaneous magnetization (| and its second derivative in the 
minimum V^g(M) (that corresponds to inverse susceptibility obtained by fitting 

either 2 or 3 histograms. 

From Table one observes that 

1) The best results for both M and V^'4(M), as far as statistical errors are concerned, 
come from the direct measurement (apparently because the region \(f\ < M is lost for 
the fitting). 

2) Fitting 2 histograms leads to results compatible with direct measurement, for all 
lattices, but with unsatisfactorily large statistical errors. 

3) Fitting 3 histograms leads to considerable reduction of statistical errors. How- 
ever, with the exception of the two largest lattices, there is a serious deviation of M 
and V^g(M) from the direct measurement. This is a manifestation of a systematic 
error associated with the finite cutoff, due to the smallness of the correlation length at 
J = J2 (corresponding to the third histogram) on our smaller lattices. 
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The improved fitting procedure 

Thus both fitting procedures in the Table |l| are far from satisfactory. So we have 
designed the following procedure that uses the J = data to significantly reduce 
statistical errors of the 2-histogram fit: one should fix the minimum of V^s at the value 
of M obtained from the direct measurement. This leads to 2-parameter fit instead of 
the 3-parameter one, with corresponding reduction of statistical errors. The results are 
collected in Table 0. 

One observes that: 1) the quality of fits, as indicated by x^? is very high: there is no 
discrepancy between Monte Carlo histograms and the fit, other than statistical noise; 2) 
V^q{M) is completely consistent with the direct measurement; 3) Fitting 5 histograms 
for our largest lattice without fixing the position of the minimum of V^s (Fig. |I| and 
the last column in Table 0) provides both M and V^g{M) completely consistent with 
the direct measurement. All parameters are consistent with the 2-histogram fit. 



CcS = ^d^(fd^ip + rLf'^ + ULp'^ + w^^ . (5) 



Extraction of universal parameters 

In the FM case the extraction of universal, i.e. dimensionless, parameters of V^s from 
the data becomes less trivial than in the PM case. 

Having measured r, w, u from the fits and the field renormalization constant Z from 
the small-momentum behavior of the propagator, we obtain the effective Lagrangian 
in the form 

1 
2 

Then we change the scale of y?, introducing the renormalized field = Z^^^^ip: 

C-eS = ^d^(fRd^yDR + ZrLp\ + Z'^ulPr + Z^w(p%. (6) 

The coefficients in front of ip'^, and cp^ have, correspondingly, dimensionalities m^, 
m and 1. In the PM case it was natural to choose y/2Zr as a scale factor, and to use it 
to render the 4-point coupling dimensionless, obtaining two dimensionless parameters: 
(74 = Z'^u/ y/2Zr and Qq = Z^w. 

In the FM case the coefficient in front of ip'ji does not determine the correlation 
length any more. Moreover, r is determined with very large statistical errors. So we 
have to find something different that could serve as a scale factor. Spontaneous mag- 
netization seems to be the best choice: (Iv^rI)^ = = Z~^M^ has the dimensionality 
of mass and very small statistical error. 

Thus the coefficient in front of (p^ provides a dimensionless parameter Z^u/M"^ 
(included in Table H). Another interesting dimensionless parameter is obtained from 
the mass: 

G - ^/M^ = (7) 

Its remarkable property is that, being a special universal combination of critical am- 
plitudes, it can be measured, to high precision, without any fitting (Fig. H): 

G = xr'M-\ (8) 
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Figure 3: Dimensionless parameter G = Z^^^^/V^M~'^ = x^~^M~^, as a function 
of inverse correlation length. The diamonds are our Monte Carlo results, the square 
is from and the solid line is from approximants given by Liu and Fisher ^ (see 
Appendix) . The line has an uncertainty of order 3% in the overall factor. If one 
combines our points with m,-dependence from Liu and Fisher's approximants, one gets 
in the scaling limit G C-{f{y^B~^ ^ 5.0 . 



lim G 



C-{f,-)-''B-'. 



(9) 



(See Appendix for notation). More dimensionless parameters are obtained from the 
third and fourth derivatives of V^s in the minimum (Table 0), such as 



K 



V"'M 
6V" ' 



(10) 



Expansion around the minimum 

Consider the situation in the vicinity of the minimum of Ves- One can write (p{x) = 
M + (p{x), where is a small deviation. Then for ip we get 

Aff = lz-%^d,^ + \v"^' + ^V'"^' + ^V""^' + . . . , (11) 
where all derivatives of Ves are taken in its minimum. After changing the scale, (p = 

^eff = ld,^Rd,^R + \ZV"^1 + j^Z^'^V'^l + \^Z'V""^\ + ..., (12) 
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Figure 4: Dimensionless parameter g^, deGned by eq. (|T^), as a function of inverse 
correlation length. Diamonds are ours, triangles are from and the square is from 



Now we obtain the mass, 



m = C^ = VZV" 
and dimensionless 3-point and 4-point couplings: 



^3 



3! m^/^ 
1 Z^V"" 



6 V 

1 Z^/'^V"" 



V" 



4! m 24 y/V" 

^(3) .^A 



(13) 

(14) 
(15) 



They are in the following correspondence with g^^ and gj^ computed by J.-K. Kim 
and A. Patrascioiu H: 
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9i 



' ~(3) 

:9r , 



24 



# + 3(^)? 



(16) 
(17) 



As can be seen from Fig. ^, our values of g^ are somewhat lower than than those 
obtained in @] by direct measurement of the 3-point correlation function. As a consis- 
tency check, we have computed g^ and g^ for two largest lattices directly from 1-particle 
irreducible 3- and 4-point functions, and obtained 



58^, 
743, 



/? = 0.2232: ^3 = 2.41(7), ^4 = 5.4(7); 
/? = 0.2227: ^3 = 2.35(8), ^4 = 5.1 ± 2.1, 



:i8) 
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in complete agreement with Table ^. 



Compact formula for V^ff 

We see that quite a few dimensionless parameters of Vcs can be constructed. However, 
one would like to have Ves written down in a concise form, so that its shape can be 
better understood. For the PM phase it could be written as 

Veff = ^m^^R + mgiip% + g^^% (19) 

with the scale parameter m and two dimensionless (and thus universal) parameters g^^ 
and 

In the broken phase one can take as a starting point the following observation. If 
the constant term in V^g is chosen in such a way that Vce = in the minimum, it has 
double zeros at lpr = ^Mr, and must be proportional to {ipj^ — M|.)^. Thus it must 
take the form 

V^si^R) = i^l - Ml)\aipl + hMl), (20) 

with two dimensionless parameters a and 6, that can be expressed via G and k (P|,pi]D, 
leading to 

V.^{Vr) = ^G'i^l - Mlf [2Ml + {2k - - Ml)] . (21) 

Unlike the symmetric phase, where we had the range of ^ sufficient to discuss extrap- 
olation to ^ oo, our results for the broken phase indicate the existence of the finite 
cutoff effects, but are clearly not sufficient to extrapolate to ^ — oo. So we can only 
write down the result for our largest ^ (the last columns of Tables |1| and 0) : 

G = 5.27(9), K = 0.998(14), (22) 
Kff = 1.73(6) ■(^^-M2)2[2M2+ 1.00(3) -(v^^-M^)]. (23) 

The hope is that these numbers are already close to the continuum limit. It is inter- 
esting to observe that Ves turns out to be proportional to {ip"^ — M^Y{ip'^ + M^). This 
sheds some light on the negative sign of ip'^, as {ip"^ — l)^(v3^ + 1) = 1 — (p'^ — ip"^ + (p^. 



Comparison with analytical results 

The only analytical results on the effective potential in the broken phase we know of 
are from the e-expansion of the scaling equation of state (see [0]), and those obtained 
recently by Berges, Tetradis and Wetterich using a method based on an exact flow 
equation for a coarse grained free energy. In their notation, 

URipR) = V,s{^r), PR = ^^l, (24) 
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Figure 5: The plots of the "Ising equation of state", J{^) = dVes{ip)/dip, with effective 
potentials computed by different means, all of them normalized in such a way that the 
minimum of V^g is at M = 1, and V^^ = 1. The curves 0-3 are from the e-expansion 



of the parametric representation of the equation of state in orders . . . e^, respec- 



tively. The curve W is from and the curve V represents our result (P^, which 
reduces in this normalization to J = — '2(p'^ + Sip^). 



they obtain 



POR 
i>R 



This translates, in our notation, to 



G 

93 



OR) 



POR 
U'HipOR) 



'^Por) 
1.08 , 

2.54 . 



61.6, 



107. 



5.55 



(25) 
(26) 

(27) 

(28) 
(29) 



These values are included in Figs. ^, ^ and are 5-8 % higher than our Monte Carlo 
results. Berges, Tetradis and Wetterich give also a complicated formula that serves as 
an approximant for J{<^) = Kfr(v^)- The corresponding curve is shown in Fig. ^, and 
goes quite close to ours. Also included are three known terms of the e-expansion of the 
parametric representation of equation of state 0. 
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Conclusions 

Here are our main conclusions about the universal properties of the effective potential 
for the broken phase of the 3D Ising model in the scaling region: 

(1) The form np"^ + mp^ + unp^ provides a very good approximation for V^Q^ip) 
(Fig. p, while without the y?^ term this is not achievable (Fig. 0) — the same result 
as in the symmetric phase fll. 

(2) Not only r, but also u turns out to be negative (Table 0). 

(3) Quantitative results can be summarized by eq. ([23|). 

(4) The combination of critical amplitudes G = C~{fi)~^B~'^ (Fig- 0) seems to be a 
quantity very suitable for a precise Monte Carlo computation, especially using improved 
estimators, as it does not suffer from uncertainties in and critical exponents. 

It should be noted that, thinking about V^fj in terms of Taylor expansion around 
the minimum, one would need at least 4 parameters to describe our data (M, V^q{M), 
V^g{M) and V^'4"(M)), and the symmetry ip —ip would be lost. Thus the possibility 
to approximate Ves with np'^ + uip^ + wip^, that has only 3 parameters and respects the 
symmetry, is by no means trivial, and demonstrates, once more, a special role of ip^ in 
three dimensions. 
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Appendix 

In Table |1| and Fig. |^ we have included, for comparison, the values obtained from 
approximants given by Liu and Fisher P|: 



x(t)=C-|t|-^-2.71, 



/i 
C 



■ 0.2502(8) 
= 0.220(4), 
1.71(2), 



2iy = 1.267, 
7 = 1.2395, 
p = 0.3305. 



(30) 




B 



We use the best available estimate of the critical coupling, Pc = 0.221655 0, rather 
than the value 0.22163 originally associated with these approximants. 
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L3 


30^ 


38^ 


46^ 


58^ 


743 


P 


0.226 


0.2246 


0.2239 


0.2232 


0.2227 


Jo 


0.0 


0.0 


0.0 


0.0 


0.0 


Jl 


0.004 


0.0022 


0.0014 


0.00078 


0.00043 


J2 


0.010 


0.0056 


0.0035 


0.0020 


0.0011 


J 










0.0017 


J 










0.0025 


AT 

'config 


3 X 360000 


3 X 360000 


6 X odUUUU 


3 X 360000 


5 X 100000 


ry—1 

^(J=0) 


0.3167(11) 


0.3166(13) 


0.3149(15) 


0.3174(14) 


0.3133(32) 


JI./T — l\,r,\\ r „ 


0.44975(17) 


0.39972(19) 


0.36743(18) 


0.32700(16) 


0.2892(3) 


M from 


0.4468 


0.3962 


0.3640 


0.3235 


0.2856 


V from 


0.0365(1) 


0.02207(8) 


0.01558(7) 


0.00958(4) 


0.00597(6) 


the propagator 












X"^ from @) 


0.0374 


0.0224 


0.0158 


0.00982 


0.00600 


e = {z-^/vfi^ 


2.946(6) 


3.788(10) 


4.496(14) 


5.756(17) 


7.24(5) 




5.30(3) 


5.22(3) 


5.23(4) 


5.12(4) 


5.27(9) 




Fitting 2 histo^ 


;rams {J = Jq, 


Jl) at \ip\ > M 


Ex' 


102 


106 


111 


92 


83 




104 


98 


108 


109 


106 


M 


0.4491(8) 


0.3993(13) 


0.3672(7) 


0.3271(7) 


0.2899(7) 


V" 


0.0353(13) 


0.0217(12) 


0.0154(5) 


0.0096(3) 


0.0063(3) 




Fitting 3 histograms (J = Jo, Ji, J2) at 93 > M 


Ex' 


188 


187 


194 


209 


109 




178 


175 


187 


188 


112 


M 


0.4470(4) 


0.3979(4) 


0.3666(9) 


0.3281(5) 


0.2894(9) 


V" 


0.0320(5) 


0.0204(6) 


0.0149(3) 


0.0102(3) 


0.0060(3) 



Table 1: Monte Carlo numerical results, computed by direct measurement of some 
observables, as well as by two different fitting procedures. Values of spontaneous 
magnetization and inverse susceptibility from approximants proposed by Liu and Fisher 
(see Appendix) are listed for comparison. A^config is the number of configurations used. 
V" = V^g(M). E x' is the sum of x' for all histograms used in the fit (it is minimized 
by the fit). Abins is the total number of bins in these histograms. The numbers in 
parentheses are standard deviations of the last decimal digits. 
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L3 


303 


383 743 


743 






Fitting 2 histograms {J = Jq, Ji) at \(p\ > M, 


Fitting 5 






keeping minimum of Ves{ip) fixed at ip = M 


histograms 








at |.^| > Af 



Ex' 


103 


106 


111 


92 


85 


200 


-^bins 


104 


98 


108 


109 


106 


204 


r 


-0.00077(21) 


-0.00084(13) 


-0.00071(8) 


-0.00036(8) 


-0.00040(10) 


-0.00038(4) 


U 


-0.0186(9) 


-0.0120(7) 


-0.0091(5) 


-0.0078(6) 


-0.0042(10) 


-0.0044(3) 


W 


0.0676(12) 


0.0612(13) 


0.0581(11) 


0.059(2) 


0.052(3) 


0.0532(7) 


M 


0.44975 


0.3997 


0.3674 


0.3270 


0.2892 


0.2893(6) 


V" 


0.03625(26) 


0.02211(15) 


0.01555(9) 


0.00955(7) 


0.00599(12) 


0.00599(13) 


m = VZV" 


0.3383(12) 


0.2643(9) 


0.2222(7) 


0.1735(7) 


0.1382(14) 


0.1382(15) 


V'" 


0.537(4) 


0.353(3) 


0.265(2) 


0.1870(23) 


0.123(3) 


0.1239(10) 


V"" 


4.48(7) 


3.23(6) 


2.60(4) 


2.09(5) 


1.47(8) 


1.497(10) 






The following parameters are dimensionless: 






-2.90(14) 


-2.37(14) 


-2.17(12) 


-2.29(19) 


-1.6(4) 


-1.72(11) 


ge = Z^w 


2.13(4) 


1.93(4) 


1.86(4) 


1.85(5) 


1.70(10) 


1.73(2) 


93 


2.55(3) 


2.43(3) 


2.39(3) 


2.41(4) 


2.27(9) 


2.29(2) 


h 


5.50(10) 


5.08(11) 


4.92(9) 


4.99(14) 


4.5(3) 


4.60(8) 




The following parameters are dimensionless and do not contain Z: 


K = V"'M/6V" 


1.110(15) 


1.065(16) 


1.045(13) 


1.067(20) 


0.99(4) 


0.998(14) 


V""M^/2AV" 


1.041(23) 


0.972(24) 


0.942(20) 


0.975(30) 


0.86(7) 


0.87(2) 




0.844(5) 


0.857(5) 


0.863(4) 


0.857(6) 


0.878(11) 


0.875(4) 



Table 2: Numerical results obtained by the fitting procedure that we find most suc- 
cessful. For the largest lattice, fitting all 5 histograms gives the same results (the last 
column), providing an additional consistency check. All derivatives of Ves are taken in 
the minimum. 
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